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Abstract 

In this paper, by using an analogy among quantum mechanics, electromagnetic 
beam optics in optical fibers, and charge particle beam dynamics, we introduce the 
concept of coherent states for charged particle beams in the framework of the Ther- 
mal Wave Model (TWM). We give a physical meaning of the Gaussian- like coherent 
structures of charged particle distribution that are both naturally and artificially 
produced in an accelerating machine in terms of the concept of coherent states 
widely used in quantum mechanics and in quantum optics. According to TWM, 
this can be done by using a Schrodinger-like equation for a complex function, the 
so-called beam wave function (BWF), whose squared modulus is proportional to the 
transverse beam density profile, where Planck's constant and the time are replaced 
by the transverse beam emittance and by the propagation coordinate, respectively. 
The evolution of the particle beam, whose initial BWF is assumed to be the sim- 
plest coherent state (ground-like state) associated with the beam, in an infinite 1-D 
quadrupole-like device with small sextupole and octupole aberrations, is analitically 
and numerically investigated. 
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1 Introduction 



It has been recently pointed out that the charged particle beam transport through a 
generic optical device can be described by means of the so-called Thermal Wave Model 
(TWM) for particle beam dynamics QTJ] . According to this approach the transverse (lon- 
gitudinal) dynamics of a particle beam is described in terms of a complex function, called 
beam wave function (BWF), whose squared modulus gives the transverse (longitudinal) 
density profile of the beam 0-||, f§ (@,@],0,@]). The BWF satisfies a Schrddinger- 
like equation where Planck's constant is replaced by the transverse (longitudinal) beam 
emittance 0-0,0 (i,@,Q,i). 

This model has been successfully applied to a number of linear and nonlinear problems 
of particle beam dynamics [l|]-[f§. In particular, it has been used in the transverse dynam- 
ics for describing the optics of a charged particle beam in a thin quadrupole-like lens with 
small sextupole and octupole deviations |J, and for estimating the luminosity in linear 
colliders |J . TWM has also been applied for describing the self-consistent interaction of a 
relativistic electron-positron beams with a collisionless, overdense plasma |2j , reproducing 
the main results for the beam filamentation threshold and the self-bunching equilibrium 

TWM seems also useful for describing longitudinal particle bunch dynamics in both con- 
ventional accelerators and plasma-based accelerators for which coherent instability, in 
terms of modulational instability, and soliton formation has been investigated 0,@,0. 
Recently, TWM has been used for describing the longitudinal dynamics of a relativistic 
charged particle bunch in a circular accelerating machine, when the radio -frequency (RF) 
potential well and the self-interaction (wake fields) are present || as well as synchrotron 
radiation damping and quantum excitation (photon noise) are taken into account M. 

The paper is organized as follows: in section 2 we briefly review the basic features 
of the Thermal Wave Model and introduce the concept of coherent state for charged 
particle beam dynamics in analogy with the light beam optics. In section 3, we extend 
the perturbative approach previously discussed for a thin lens 0], M to investigate the 
stationary propagation of a charged particle beam through an infinite 1-D quadrupole 
with small sextupole and octupole deviations. The initial BWF is assumed to be in the 
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ground-like state of the harmonic oscillator (pure Gaussian profile). Numerical results for 
the BWF obtained with this approach are presented and discussed. Finally, in section 4 
we give our remarks and conclusions. 



2 The concept of coherent state in particle beams 

Let us consider a relativistic charged particle beam travelling along z-axis with velocity 
(3c (/3 ~ 1). Denoting by x the transverse coordinate of a single particle in the beam and 
by p its conjugate momentum, we consider the following single-particle hamiltonian 

2 

H = P - + U{x,z) . (1) 

where H has been made dimensionless dividing by the quantity m , yf3 2 c 2 (m and 7 being 
the particle rest mass and the relativistic factor, respectively). In ([I]) U(x,z) stands for 
a dimensionless potential filled by the particle. 

2.1 Main features of TWM 

As stated before in the Thermal Wave Model the collective behaviour of a charged 

particle beam is ruled by a Schrodinger-like equation, which can be obtained from the 
single-particle description by means of the following Thermal Quantization Rules (TQR) 

M, I, 

d - d 

x — ► x , p — > p = —ie— , and H — > H = ie— . (2) 

ox oz 

where here e stands for the transverse beam emittance. Consequently, prescriptions (0) 
give immediately the following evolution equation for the BWF 

d e 2 d 2 

ie—y{x,z) = -- — V{x,z) + U{x,z)y{x,z) . (3) 

Provided that the normalization condition \^{x, z)\ 2 dx = 1 is satisfied, the trans- 
verse number density of the beam particles A(x, z) is given by 

A(x,z) = N\V(x,z)\ 2 , (4) 

where N is the total number of particles. Eq. (^) is formally identical to the paraxial wave 
equation which describes the optical beam propagation in inhomogeneous media with 
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refractive index profile given by U(x, z). In this analogy the inverse of the wave number 
is replaced by the beam emittance. 

It has been shown in JI[ that in the simplest case of a relativistic charged particle 
beam crossing a pure quadrupole (aberrationless lens) of focusing strength k\ > , the 
BWF satisfies the following parabolic equation 



■ 9 t, ^ 



e 2 d 2 1 
2*(z, z) + -h x 2 ^(x, z) 



2 dx 2 ~ v "' '' ' 2 

It is easy to show that Eq. (|5]) admits the following orthonormal discrete modes 



(5) 



x exp 



x 



„2 



1 



x 



~ [2tt a 2 (z) 2 2n (n!) 2 ] 1/4 
i(l + 2n)(f)(z 



H„ 



x 



y/2, 



a(z 



with n = 0,1,2, .. 



4a 2 (z) ' 2ep(z] 

where H n are the Hermite polynomials, a(z) obeys to the following envelope equation 



(6) 



a" + ki o + 



4a 2 







(7) 



and 



1 

P o 



Aa 2 



(9) 



where each prime denotes the derivative with respect to z. Within the similarity with 
electromagnetic beams, ^ n (x,z) play the role analogous to the one played by Hermite- 
Gauss electromagnetic modes, a(z) describes the beam caustic (i.e. beam envelope), and 
p{z) represent the wavefront curvature radius. By introducing the r.m.s. cr x (z) associated 
to a general BWF satisfying (|3]) as (the mean value of x is assumed equal to zero) 

x) l 1/2 



a x (z) = (x 2 ) 1/2 



x \ty(x,z)\ dx 



(10) 



(the quantum-like expectation value), it is easy to see that o~(z), appearing in (Q)-( 
coincides with the r.m.s. of the fundamental mode ^o(x,z) 

1/2 



o [z 



x \^fn(x,z)\ dx 



In addition, we can also define the expectation value for the transverse linear momentum 
associated to *f? (x,z) 

1/2 

a p {z) = e I — r/.r . (12) 
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By following particle accelerator physics language, cr(z), p(z) and <r p (z) can be expressed 
in terms of some optical parameters, called Twiss parameters a(z), (3(z), and 7(2) || 

a 2 (z) = e(3(z) , (13) 

— = (14) 
a 2 p {z) = e 1 {z) . (15) 

It is suitable to introduce the following matrix 

- ( t) M ) ' (16) 

whose determinant, as it is well-known ||, is constant with respect to z, namely 

d [detf (z)] /dz = . (17) 

Consequently, we choose for our convenience, without loss of generality, 7/? — a 2 = 1/4. 
Consequently from (pT3|)- (p~5|) follows that the determinant of the matrix 

~ a{ % a [ {Z) ) • (18) 
v ' \ —a(z) a (z) a (z) / v ' 



is an invariant, namely 

2 2 / /\2 

rr rr — inn ) 

; 

It is easy to prove that 



o-;cr 2 - (aa'Y = — = const. . (19) 



a(z)a (z) = / * (a;,z) o Vo(x,z) dx = ( ) . (20) 



Consequently, from flI3p follows that 



(x 2 ) (f) ~ (^^) 2 = € ~ , (21) 
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which is formally identical to Robertson-Schrodinger uncertainty relation [[ITJ. [I Ij] for 
partial cases when (x) = (p) = 0. Furthermore, within the framework of TWM, (EH]) is 



a quantum-like version of the well-known Courant-Snyder invariant |12[]. Note that (19[) , 
or equivalently fl21~D , gives immediately the usual form of the Heisenberg-like uncertainty 
principle [0] which is analogous to Heisenberg uncertainty relation in quantum mechanics 
[L3| (again for (x) = (p) = 0) 



a P a> € - . (22) 



By introducing, for the general hamiltonian (|1|), the total averaged energy associated to 
the transverse motion of the particles 



V*(x,z)HV(x,z) dx = - / dx+ / Z7 |*(ar,«) 

(23) 



I cfx 



the following virial equation 



cLs 2 dx 
and the following energy-variation equation 

d£ r+^dU |2 



d2 ^ ) -4f-2<xS , (24) 



.. ,*(x,z)r cfe (25) 

(12 J-oo (72: 

hold. For a quadrupole-like potential (see Eq.(|5|)), Eqs. (p4|) and ( p5| ) show that £ is a 
constant of motion (dS/dz = 0), and , in particular, for \1/ = \l/o recover the envelope 
equation (|7|). The equilibrium solution of (0) (d 2 a(z)/dz 2 = 0), namely 

ff -27ir • (26) 

plays an interesting role. In fact, (^) prescribes that if we prepare the system, before 
entering the quadrupole-like device, with given a , e and k\ in such a way to satisfy 
(p6|), the evolution of the beam through this device is performed with the transverse size 
fixed at the equilibrium value (<j(z) = ctq). In other words, ( p6f ) corresponds to an initial 
beam configuration with the Twiss parameter a(z) = 0, i.e. zero-divergence of the beam, 
or equivalently to a BWF with a wavefront whose curvature radius is infinity. Thus, 
the beam divergence (the curvature radius of the BWF wavefront) is also zero (infinity) 
during the beam evolution. 



Note that for a quadrupole-like potential and for the equilibrium solution fl26|), the set of 



Hermite-Gauss modes @ reduces to the hamiltonian eigenstates of the harmonic oscillator 

') = j^pFF exp (~5 + ' (1 + 2n)Mz> ) fe) ' (27) 

where n — 0, 1, 2, 

00^) = -^ , (28) 
and the energy values given by (|23|) under the substitution of these eigenstates, coin- 
cide with the hamiltonian eigenvalues of the harmonic oscillator 

£°=(n+£)e^ . (29) 



In particular, for n = (j2TD, ( p8|) and (|29| ) give the ground-like state 

*) = [ 27r jg]V4 exp ( _ ^f + , (30) 

which is a pure real Gaussian and the lowest energy reachable by the beam is Sq = 
{l/2)ey/k[. The means (x) and (p) are equal to zero at this state of beam. In these 
conditions the uncertainty relation is minimized as 

O"0 CTpQ = - . (31) 

Eq. (^l|) holds also during the evolution of the beam, because, in addition to (^6), we 
have u p {z) = <j p o = const.. In summary, we conclude that if we initially prepare the 
beam according to the matching conditions (p6|), its evolution is ruled by a quantum-like 
behaviour in terms of BWF-ground-like state which minimizes the uncertainty relation 
and corresponds to the lowest accessible beam energy (l/2)e-y/^i- 

By introducing the connection between emittance and transverse beam temperature 



T, showed both in the conventional beam physics 0, and in TWM 0, e 2 /Aal = 
/csT/(mo7/3 2 c 2 ) (ks being the Boltzmann constant), we easily obtain 

_ k B T 

where ( PB| ) has been used. Consequently, the beam energy associated with the equilib- 
rium solution corresponds to the thermal energy and this result is in agreement with the 
equipartition theorem. 

As we will show in subsection 2.3, BWF ([30]) belongs to the infinite series of coherent 
state functions, labeled by a complex number a = a\ + ia2, and widely used in quan- 
tum mechanics and quantum optics fT5|-|]l7|. To further develop this point, in the next 



subsection we present the general formalism of coherent states and their main properties. 
2.2 Formalism of coherent states 



In fiber optics the complex ray formalism |L8 |, ||19f has been used to describe Gaussian 



wavepackets propagating along the fiber in frame of Fock-Leontovich parabolic equation 



20H which is a Schrodinger-like equation. In Refs. [^T[] , [22] it was proved that the complex 



rays are just the coherent states describing the electromagnetic beams in fibers. Within 
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Dirac's bra and ket formalism, the coherent state \a) is the state which represents the 
following series 

/ \ a \ 2 \ » a n , . 

|a> = exp — ) 2^n=A n ) > ( 33 ) 

V 2 J n=oVn\ 

where the number state \n) corresponds to discrete modes of the electromagnetic beam 
and a = ot\ + % a<i is a complex number. Since the following orthonormal condition holds 

(n\m) = 5 nm , (34) 

we have the normalization condition for coherent states: 

(a\a) = 1 . (35) 

A coherent state may be obtained using the unitary shift operator 

V a = exp (a a) — a* aj , (36) 

where the photon creation and annihilation operators a} and a satisfy the boson commu- 
tation relation 

1 • (37) 



a, a* 



Thus, if one acts with this operator on the fundamental mode state |0), such that 

a |0) = , (38) 

we have the normalized coherent state 

P a |0) = |a) . (39) 

Note that V a has the following shifting properties 

V ] a aV a = a + a , (40) 
Vla^V a = a* + a* . (41) 

(42) 

From these formulas it follows, that the coherent state is the eigenstate of the annihilation 
operator 

a \a) = a\a) . (43) 
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If one calculates in the coherent state the means of quadrature components Q and P 
one has 

(a\Q\a) — -^={a\a + at\a) — a/2o:i , (44) 

(a|P|a) = —j=(a\a — a)\a) = V2a 2 ■ (45) 
V2 

The shift operator has the multiplication property 

V a V p = V a+(i exp (ilm[a(3*\) , (46) 

where Im[...} stands for the imaginary part. The distribution function of the electromag- 
netic modes \n) in the coherent state 

P(n) = \(n\a)\ 2 = exp H^ 2 ) \a\ 2n (47) 
is the Poisson distribution for which the expectation value of mode number 

oo 

(n) = J^n P(n) = \a\ 2 , (48) 

n=0 

and the dispersion of the mode number 

On = (n 2 ) - (n) 2 = \a\ 2 (49) 

is equal to the same value. For coherent states the dispersion of photon quadrature 
components 

{a \( a +f)'\a) -((a|^|a)) 2 = ± , (50) 



^2 J ' ' 2 ' 

Mp-yflV) -(<«|-#l«))' 2 = \ , (51) 
are equal and there is no correlation of these quadrature components 

(a\ | a + a\ a — a^j \a) — 2(a\a + a)\a) (a\a — a)\a) = , (52) 

where { , } denotes the anticommutator of two operators. In the coherent state repre- 
sentation the wave function of discrete mode state has the simple form 

... a n ( \a\ 2 \ . . 

(«F> = -r=fexp — . (53) 



As one sees from the decomposition Q33|), the state exp(|a| 2 /2)|a) is the generating func- 
tion for the number state \n). 

According to ([50]) and ([H]) the products of the dispersions satisfies the condition of mini- 
mizing the Heisenberg uncertainty relation. The evolution of the coherent state \a) along 
the beam path due to the hamiltonian of the oscillator H = u (a^a + 1/2) produces again 
the coherent state with extra phase factor 



cxp 



— iusz [ a) a + - 



| a) = exp 



.toz 



| a exp (— iujz)) . (54) 



From that point of view the coherent state is considered as the closest to the classical state 
of the oscillator since it represents the wavepackets whose center moves along the classical 
trajectory in the phase space and the width of the packet is minimal and consistent with 



the uncertainty relation [|15]-[17|, []23|, |24 



In the next subsection we will discuss the properties of the coherent states in the coor- 
dinate representation in order to describe the coherent charged particle beam transverse 
dynamics within the framework of TWM. 

2.3 Physical meaning of coherent states for charged particle 
beams 

One of the needs that we understand to introduce the notion of coherent states for charged 
particle beams is connected with the experimental possibility of producing, by acting with 
an external electromagnetic force, a shifting off-axis of the beam propagation coordinate. 
For example, in an accelerating machine this can be done by means of some devices, 
such as kickers, RF cavities, etc. [[25]-[27|, which are able to produce a displacement 



of the transverse beam distribution center from its stationary position (design orbit or 
synchronous particle orbit) to another neighbouring path. This means that the center 
of transverse beam space-distribution has a shift and this effect allows us to naturally 
introduce, within the framework of TWM, the coherent state representation. In fact, let 
us suppose that the beam is in the ground-like state ^/q(x,z) given by fl27|) for n = 0, 
which describes the equilibrium state of the particle beam travelling along the stationary 
orbit |[28||. In order to take into account the displacement of the transverse distribution 
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center we introduce the following modified potential well 

1 



U(x,z) 



-k\x 



F x {z)x + F 2 (z) , 



(55) 



where F\(z) and F%{z) account for the external electromagnetic forces. Since U(x,z) is 
still quadratic with respect to x, from general properties of Schrodinger equation, the 
initial ground-like state \I/q continues to be Gaussian, but its center (wavepacket center) 
is now shifted off-axis. In order to explicitely find this transformed BWF, let us consider 
the following transformations on @ with U given by (|55|) : 



ty(x, z) = exp 



i 

e Jo 



Fo(z') dz' 



e(x,z) 



(56) 



which allows us to get 

ie^-Q(x,z) 



e 2 d 2 /l 

-— Q(x, z) + i^-h x 2 - F x {z) x ) 9(x, z) 



(note that |\I/(x, z) 
Gaussian solution 

®q{x,z) 



\®(x,z)\). It is very easy to prove that 



(57) 

has the following 



1 


1/4 




exp 


2 7T a 2 _ 





(x - x (z)) I , . , . 

2 1" -P0{Z) X ~ WOW 

4 a e 



(58) 



where <7o still satisfies the equilibrium condition (|26[ ) and the functions xo(z), po(z) } and 
Sq(z) satisfy the following differential equations 



x' ' + hx = Fi(z) 
Po + hPo = F[{z) 



and 



60 = Ye 



ki--9- 4- 

1 A 9 

4^o 



(59) 
(60) 

(61) 



Some considerations are in order: 

i) . Eqs. fl5"9"D and (^) coincide with the constraints for the following canonical transfor- 
mation x,p — > x,p, with 

x = x — xo(-2) (62) 
P = p-Po(z) (63) 

which transforms the hamiltonian H = p 2 + x 2 — -^1(2;) x into H = p 2 + x 2 . 

ii) . Consequently, Xq andpo represent the space-coordinate and the momentum-coordinate 



11 



shift, respectively. For the following it is suitable to introduce the following complex 
dimensionless shift 

Xo(z) ,cr p (z) 



am 



2a 



+ i- 



aAz) + za 2 m 



(64) 



iii). Eq.(^) (Eq.(^)) can be thought as classical motion equation of unit mass particle 
which feels the restoring force —k\Xo (— &ipo) and the external force F\(z) (F{(z)). 



iv). It is very easy to see that, by virtue of (p6|), solution (58) still minimizes the uncer- 
tainty relation, and, consequently, describes a coherent state associated to the shifts x 
and Pq. 



v). According to the definition (^), solution (58) gives the following coherent-state- 
particle-distribution Aq(x,z) 



A (x,z) = N \e Q (x,z) 



N 



: exp 



(x - 2a a 1 (z)Y 
2al 



(65) 



Thus, 
where 



can be cast in the following form 



2 7T a$ 



1/4 



exp 



—x 

ToJ 



exp 



a(z) x \ a (z)\ 2 a 2 (z) 
~^~ 2 2~ 



z = 



—oc 1 {z)F 1 (z) — 



dz 



exp [i9(z)\ , 
(66) 

(67) 



Note that, according to the formalism used in subsection 2.2, the BWF Qq(x,z) is a 
coherent state produced by the complex shift a{z). In the Dirac's formalism, let us both 
denote this coherent state as |a), and introduce the notation 



(x\a) = 6o(x, z) 



For the sake of simplicity we will consider now these states at z = 0, namely 



(68) 



(x\a ) = Qq(x,0) = $q, (x) 



1/4 



exp 



— x 
4^ 



exp 



OtQ X |CKq| 



(69) 



where a = a(0) = ai(0) + za 2 (0) = a w + ia 2 o- These functions have the following 
properties pT5[|- [|T7[| 

Kl 2 |A>| 2 



/+oo / 
®*a ( x )®Po( x ) dx = exp - 
-oo \ 



(70) 
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which means nonorthogonality, and 



] r+oo r+oo 

da 10 da 2 o \a ) (a \ = 1 , (71) 



oo J —oo 



71 

which implies the completeness condition, i.e. any state |$) with the wave function 
may be represented as a superposition of the coherent states: 



I r+oo r+oo 

|$) = - / da w / da 20 (a |$) |a ) , (72) 

7T J —oo J —oo 

where the BWF in coherent state representation (ao|^) is given by the overlap integral 

/+oo 
K (xMx) dx . (73) 
-oo 

Since the overlap integral of two different coherent states is not equal to zero these states 
form an overcomplete set of functions. The particle beam coherent state $ Q0 (x) is a 
normalized eigenstate of the annihilation operator 

^Tffe + ^S ' (74) 

i.e. 

a § ao (x) = a ® ao (x) . (75) 

Since TWM is described for small intensities by the linear Schrodinger-like equation, 
any BWF associated to a charged particle beam can be represented as a continuous 
superposition of coherent states, or discrete superposition of Fock mode states. 

3 Charged particle beam dynamics through an infi- 
nite 1-D quadrupole with small sextupole and oc- 
tupole deviations 

In this section we describe the transverse dynamics of the charge particle beam while it is 
travelling through an infinite 1-D quadrupole with small sextupole and octupole deviations 
(aberrations). We consider here the same kinematic and geometrical assumptions made 
in the introduction of section 2. 

Without aberrations the beam is assumed to satisfy the matching condition (|26|) (un- 
perturbed beam). Thus, according to the results of section 2, if the initial (z = 0) 
transverse beam distribution is purely Gaussian, i.e. 
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the unperturbed beam evolution at any z > is described by a particular coherent state 
which is the ground-like state fl30|) that we can obtain from (|58|) for the special case 
x o — Po — or, equivalently, from for a = : 

= [2^ eXP (^) eXP ("^) • (7?) 

In order to take into account small sextupole and octupole aberrations, we perturb the 
(Bp into the following equation: 

d e 2 d 2 1 

ie— #(z,z) = —2Q-2 ^(x,z) + -k 1 x 2 ^(x,z) + V ^(x,z) , (78) 

with the initial condition (f76|) , where the perturbation 

V-( x ) = I A; 2 x 3 + i A; 3 x 4 , (79) 

accounts for the aberrations. In particular ^ k 2 x 3 k 3 x 4 ) corresponds to the 1-D 
sextupole (octupole) potential term. Note that V(x) given by (ff9|) can be considered a 
small perturbation provided that the following conditions 

k 2 (r /3k l « 1 and k 3 a 2 /12h « 1 (80) 

hold. In the next subsection we will give an approximate solution of the problem ( [78] ) 
with the specifications (f76|) , (ff9|) , and (jS0|). 



3.1 Analytical results: perturbative approach 

By using the standard perturbation approach we write the following expansion of the 
ty(x, z) in terms of the eigenstates ty° m 

+oo 

^(X,z)=j2cm(z)^° m (x,z) . (81) 
m=0 



The substitution of Eq. ( plD into (|78f) yields the infinite set of equations 

dc + °° 

ie-^ = E c m{z) (n\V\m) , (82) 

m=0 



where |n) = To solve the set of equation (821) we consider the first order correction to 



the BWF in the case in which c m (0) = S mt o, according to the initial condition (|76|) . If we 
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write c m (z) = 5 mt0 + c^(z) (conditions ( p0| ) here correspond to the condition \c l m \ « \) 
Eqs. (pi) give 



■ d 1/ X 



(nl^lO) 



where 



(n|V|0)= / **(ar, i 2)V r (a;)*g(ar,«) 

J — oo 

Hence, the coefficients c\(z) are given by 



c}(z) 



Co(Z 



cl(z) 



c](z) 



■ 3 t 



2^/2 



v exp 



l- 
2 



■ (* 
sin - 

\2. 



-i-fi exp sin(£) 



(83) 
(84) 

(85) 
(86) 
(87) 



12^2 
1 



f exp 



.3e 



sm 



'3^ 



i9) 



where we have introduced the dimensionless parameters v = fc20"o/3fci and \i = &30o/12A;i, 
and the dimensionless length £ = \fk[z. The non-normalized perturbed BWF is then 
obtained by the superposition ([H]) for four terms only: 



[1 + <%(z))*° (x, z) + cj(z)*?(x, z) + cl(z)y° 2 (x, z) 



+ c\{z)^{x,z)+c\{z)^l{x,z) 



(90) 



Note that now the BWF ^(x, z) given by ( |90|) is not Gaussian anymore and, consequently, 
does not represent a coherent state anymore. In fact, the aberrations introduce a defo- 
cusing of the particles which produces a distortion of the particle beam distribution with 
respect the unperturbed Gaussian profile. 

In the next subsection we analyze more accurately this distortion on the basis of 
numerical results. 

3.2 Numerical results 

According to the perturbation theory results obtained in the previous subsection, the 
parameters /1 and v represent a measure of the distortion due to sextupole and octupole, 
respectively, once the multipole distortion is defined as the ratio of the sextupole (oc- 
tupole) aberration strength to the quadrupole strength. Provided that fi « 1 and 



15 



v « 1, the distribution in configuration space is proportional to \^\ 2 , where \1/ is given 
by (0). In Fig.l we have plotted the normalized transverse density profile versus the di- 
mensionless transverse coordinate x/^/2a for increasing values of £ and for /i = .005 and 
v = .05 (see Figs.la — lf). The dashed lines give the starting distribution. We can clearly 
see a weak distortion of the particle distribution as the beam propagates throughout the 
optical device, which is quite evident in Figs. lb, Id, and 1/: the solid lines represent the 
distorted distributions for £ = 15.5, 21.5, 27.5, respectively. As £ increases monotonically, 
the distortion increases and decreases alternatively. In particular in Figs.la, lc, and le 
(£ = 12.5,18.5,24.5, respectively), the distortion is negligeable (negligible distortion), 
while Figs. lb, Id, and If show significative distortion. This effect is due to the fact that, 
according to (0), the BWF is given by a superposition of five modes which interfere each 
other. To clarify this point in Fig.2 we have plotted the evolution of the normalized 
distortion of the transverse density profile as a function of the dimensionless longitudinal 
coordinate £. The normalized distortion is defined here by Aa/a , where Act = cr(£) — cr 
is the deviation of the effective transverse beam size (r.m.s.) evaluated at a generic £ 
from the initial one <7q. In Fig.2 we have fixed v = .05 and plotted the distortion for 
/i = 0, .0015, .002. When \i = 0, the octupole gives no contribution and the BWF is a 
superposition of three modes only (see Eq.(pU|)). From Fig.2 is quite evident that the 
oscillating behaviour of the distortion increases if we increase \i. 

In Fig.?> we have plotted the evolution of the distortion fixing /i = .005 and for v = 0, .05. 
When v = 0, we have no contribution from the sextupole, but, when both sextupole 
and octupole deviation are taken into account, we clearly see the enhancement of the 
distortion of the beam profile. 

4 Remarks and conclusions 

In this paper we have introduced the concept of coherent state for charged particle beams 
within the framework of TWM 0-0. 

We have shown that the coherent Gaussian-like structures for charged particle beams, 



naturally or artificially generated in particle accelerators |25|- ||271| , correspond to examples 
of beam wave functions which have a formal definition identical to those ones are widely 
used in the quantum mechanics and quantum optics [15]. 
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We have used the analogy among quantum mechanics, electromagnetic beam optics in 
optical fibers, and particle beam dynamics in the framework of TWM in order to introduce 
the concept and to understand the physical meaning of coherent state for a charged particle 
beam. We have shown that such a kind of state is associated with the complex parameter 
a, introduced in section 2, which describes the shift of the beam distribution center from 
the optical axis. In particular, we have proved that when the charged particles of a purely 
Gaussian beam, travelling in an accelerating machine, are shifted from their path by the 
action of some external transverse electromagnetic forces H25|-|27||, the state after this 
action is just the coherent state |a) (the real part of a corresponds to the shift of the 
space coordinate x and the imaginary part of a corresponds to the shift of the conjugate 
momentum p) . 

Physically, due to the action of an external electromagnetic force, the distribution 
center results to be shifted. This shift, in turn, transforms the initial BWF modulus 
|\&g(x, z)\ into the following: \^{x, z)\ = \^q(x — Xq, z)\. In addition, since the shifting does 
not change the equilibrium condition (|26|), because it is an intrinsic relationship between 
the quantities <To, e, and k\, and since the BWF profile is still Gaussian, the uncertainty 
relation initially fixed at its minimum value e is preserved during this shift. But this 
picture, in terms of formalism and physics of charged particle beams, fully corresponds to 
the coherent state description given in subsection 2.2, and the transformed BWF describes 
a new coherent state with respect to the simplest one which corresponds to the ground-like 
state. From the experimental point of view the shift operator T> a (|36]), which produces 
these new coherent states for charged particle beam, in an accelerator can be realized by 



means of some devices such as kickers, RF cavities, etc. 251-127 



Furthermore, by taking the ground-like-state as the initial condition, we have studied 
the particle beam evolution through a infinite 1-D quadrupole-like in the presence of small 
sextupole and octupole deviations. We have shown that, in the framework of TWM, this 
evolution is governed by a 1-D Schrodinger-like equation for the BWF which have been 
solved by means of the standard time-dependent perturbation techniques. To the first- 
order, this perturbative treatment has shown that, during the evolution, the particle beam 
profile does not correspond to a coherent state anymore. This is due to the aberrations 
which produce a distortion of the initial Gaussian profile: during the beam evolution 
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through the optical device the BWF is a superposition of four Hermite-Gauss eigenstates 
only. The interference of these four states produces, during the evolution, small oscillating 
distortions of the particle beam profile around the ground-like state. Numerical estimates 
for this effect have also been given. 

In a forthcoming paper we will consider the present perturbative treatment by taking, 
instead of the ground-like state (a = 0), a general coherent state (a ^ 0) as initial 
condition. In addition, since in quantum optics there were introduced other quantum 
states which are generalization of coherent ones, such as squeezed and correlated states 
29R- [j3^1 , we will discuss their analogs in particle beam physics. 
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Figure captions 



Fig.l Normalized transverse profile of the beam density \Z2TTa \^\ 2 vs. x/\/2cr for in- 
creasing value of £ and for v = .05 and /x = .005. 
Dashed line = starting distribution, solid line = distorted distribution 
(a) £ = 12.5; (b) £ = 15.5; (c) £ = 18.5; (d) £ = 21.5; (e) £ = 24.5; (f) £ = 27.5; 

Fig. 2 Normalized distortion Aa/a vs. £ for fixed v = .05 and \i = 0., .0015, .002. 

Fig. 3 Normalized distortion Aa/a vs. £ for fixed pL = .005 and v — 0., .05. 
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